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' The dynamics in a nonlinear Schrodinger chain in an homogeneous electric field is studied. We 

^S) I show that discrete translational invariant integrability-breaking terms can freeze the Bloch nonlinear 

oscillations and introduce new faster frequencies in their dynamics. These phenomena are studied by 
direct numerical integration and through an adiabatic approximation. The adiabatic approximation 
allows a description in terms of an efi^ective potential that greatly clarifies the phenomena. 
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I. INTRODUCTION 
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. ' The study of the generalized discrete nonlinear Schrodinger (GDNLS) equation, introduced by Salerno [1| as a model 
. ^ ' providing one-parametric transition between discrete nonlinear Schrodinger (DNLS) equation and exactly integrable 
r— H ■ Ablowitz-Ladik (AL) model 0, is attracting increasing interest, due to the relevance of lattice dynamics in various 
I " i ' fields of physics, as condensed matter, fiber optics physics, molecular biology (see e.g. 0| and references therein) and 

recently Bose-Einstein condensate (see e.g Q and references therein). 
C^l One of the most interesting phenomena, which can be observed in the different discretizations of the nonlinear 

^ : Schrodinger equation is the so-called Bloch oscillations, that appear when a linear force is applied to a solitary wave 
solution. Such oscillations in the integrable AL model with a linear force have been discovered in numerically 
. found in the DNLS 0, and interpreted as Bloch oscillations, using the analogy with the solid state physics, in Q. 
• Later on, Bloch oscillations were studied in the GDNLS equation Q and in the presence of impurities 0. Bloch 
oscillations have been observed experimentally in an array of waveguides |3 , found also to exist in the case of a dark 
soliton and in a totally discretized (i.e. discrete-time discrete space) nonlinear Schrodinger equation 

In the present article we report some peculiarities of the Bloch oscillations of bright solitons in the GDNLS equation. 
. More specifically we show that the model preserves the Bloch oscillations, the amplitude of which however displays in 
' ^ ' certain cases an abrupt change when the parameter governing transition between AL and DNLS models is changed. 
(~| , Section ^] presents the model, i.e., the GDNLS equation with discrete translational invariant integrability-breaking 
terms in an external homogenous electric field, and the soliton initial conditions considered for this model. In Section 
nil we present the adiabatic approximation, that will prove to be very useful to understand the dynamics. It allows 
a description in terms of an effective potential. In Section llVI we describe the main ingredients of the dynamics: the 
?-j Bloch oscillations and the freezing of the Bloch oscillations, and provide an explanation of this features in terms of 
the adiabatic effective potential. Finally, in the conclusions (Section 0), we summarize the main results and discuss 
their relevance. 

II. THE MODEL 

The equation of motion for the system we are dealing with reads 

at 

where v is the integrability-breaking parameter providing one-parametric transition between the AL model 2] [v — 0) 
and the DNLS model {v — 1), and % is a parameter defining the strength of the linear force. In particular, at = 
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Eq. is integrable and has the exact bright sohton solution 
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where w and are constant parameters of the solutions (^o can be interpreted as the initial position of the soliton 
center, and l/w as the soliton width). 

Our purpose is to study how the u-ievva changes the soliton evolution; i.e., we consider the evolution for ^ of 
initial conditions of the form 

, sinh(2u>) 

= coshp^-^o)] 



III. ADIABATIC APPROXIMATION 



Let us start with the analysis of the problem when it is close to integrable, i.e., when \v\ <^ 1. We employ the 
perturbation theory (in the presence of a linear force it was developed in 0) and limit ourselves to the adiabatic 
approximation. This means that the term v^ipn+i + V'm-i ~ '^ipn)\4'n\'^ is considered as a perturbation of the AL model. 
As we have previously seen, for v — there are exact solutions of the form 

sinh(2w) e^2ie(Ti-?)+iip ^-i{2n+i)xt 

cosh[2w(n — 4)J 

In the adiabatic approximation we compute the time evolution of the parameters w, ^, and ip; while keeping the 
functional form of Eq. Q fixed. The equation of motion Q can be written as 

+ (1 + |V'n|^)('0n+l - i'n-l) - ^Xmpn = l^ilpn+l + ^n~l - 2ll}„)\'4}n\'^ . (5) 

with the z/-term considered as a perturbation. This approach gives the following evolution equations for the parameters 

^ = A{w,0coB[2{xt + e)]+B{w,0 (6) 
^ = C{w,S.) sin[2(xt + 0)] (7) 
^ = D{w,0 sin[2{xt + e)] (8) 



where 
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C(w,C) = ^ (1 + - sinh^(2w 
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From these expressions one can see that the coefficients A, B, C and D are periodic in ^ with period 1. The dynamics 
of 9, and w is decoupled from the evolution of the phase (p, and ip is merely slaved to the dynamics of the previous 
parameters. We have seen numerically that for a wide range of parameters A{w, ^ = and D{w, ^) = (for ^ integer 
of halfinteger it is easy to prove it analytically using that the terms of the series are odd- functions in n), therefore the 
adiabatic equations of motion reduce to 
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The corrections in B and C due to the iz-teim, B/v and [C + smh{2w)/w]/v [Eqs. pO|l - llHl ]. are plotted in Fig. 
for various values of w. One can see that B/v is (^-dependent and takes greater values for narrower solitons. On the 
other hand, the corrections in C coming from the v-teirtY are greater for wider solitons, and are nearly ^ independent. 




FIG. 1: (a) B/u and (b) [C + sm\i{2w) / w] / v as a function of ^ for various values of w. 

In the previous equations we can make the change of variable a = + after this change it results a system of 
autonomous differential equations for a and ^ 



da 
'dt 
di 
dt 



X + B{0 (16) 
C{£) sin(2Q) . (17) 



The differential equation for the orbits is 



da ^ X + B{0 

d^ C(e)sin(2a) ' ^ ' 



that gives the equation for the orbits a = a{£). Deriving with respect to time Eq. fTTjl and substituting the equation 
for the orbits, a = a{i/), we obtain a second order evolution equation for ^ in terms of an effective force, Fgff, that 
only depends on ^ 



FeffiO - ^C{Osm'[2am+2CiO[x + Bmco8[2am. (20) 
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This evolution equation (multiplying by d^/dt and integrating over time) leads to the energy conservation law 

Veff{0 = - f^dC'F^ffia , (22) 

where Vef / is the effective potential. Thus, the soliton in the adiabatic approximation can be seen as a particle moving 
in an effective potential. This interpretation is very useful to understand the Bloch oscillations and its freezing. 



IV. DYNAMICS IN THE PRESENCE OF THE //-TERM 

The evolution of soliton initial states of the form 



cosh[2 w[n — ^o)J 



is characterized by two main phenomena: the Bloch oscillations and its freezing. 



A. Bloch oscillations 




FIG. 2: (a) Bloch oscillations of a soliton [full evolution, Eq. 10] in a Schrodinger chain with v — —0.01. (b) Dynamics of the 
center of the soliton [adiabatic approximation, Eqs. Q-JSJ] for u — —0.01 (full line) and for i/ — (dashed line). Both figures 
with initial conditions wo — 1.5, = and 6o — 0, in external field x ~ —0.3. 



The dynamical localization can be understood in the adiabatic approximation (see Eq. I22|l . This approximation 
predicts that the soliton center position, ^, feels a trapping effective potential (see Fig. and the turning points 
of the motion are for 2{xt + 0) multiple of tt (see Eq. 114(1 . The adiabatic approximation gives good results for 
the amplitude and the period of the Bloch oscillations, and we also observe that the soliton width (1/w) remains 
approximately constant as the adiabatic approximation predicts. 

The v term changes the period and amplitude of the Bloch oscillations, and introduces new (fast) frequencies in the 
dynamics [Fig.E]. These differences can be seen as a consequence of the changes in the adiabatic effective potential, 
see Fig. O For i/ > the lattice sites have an effective attractive interaction for the soliton center; thus the effective 
potential displays local minima near the lattice sites. Conversely for v < the effective interaction is repulsive, and 
the local minima are located in the intersite regions. 

The evolution of the center of the soliton is correctly predicted by the adiabatic approximation (Eqs. ((13|l - H15(l ) at 
early times. Later the effects of the radiative decay and the change of shape of the soliton became important, and the 
description of these effects falls beyond the scope of the adiabatic approximation. These nonadiabatic effects become 
more important for increasing j^j. Another phenomenon that cannot be studied within the adiabatic approximation 
is the emergence of chaotic behavior [l3j | , the adiabatic approximation do not account of this effect because it reduces 
the computation of the evolution to the integration of the conservation law in Eq. (|21|l . 
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FIG. 3: Effective adiabatic potential for a soliton with initial conditions wq — 1.5, = and = in a Schrodinger chain 
with an external field x ~ —0.3 for various values of u. 



B. Freezing of the Bloch oscillations 

The main new phenomenon reported in this article is the freezing of the Bloch oscillations due to an homogeneous 
integrability breaking term (the I'-term). This freezing is a trapping of the soliton in a lattice site (or in a intersite 
region), that happens in a region of the parameter space [Fig. This is an intrinsic localization effect due to 
characteristics of the lattice chain, not a localization induced by lattice inhomogeneities . 




FIG. 4: (a) Narrow soliton (wo = 1.5 and = 10); and wide soliton {wo — 1.3 and = 0) evolutions in a Schrodinger chain 
with u = 0.02 in an external field x = 0.3, both with 9o = 0. (b) Effective adiabatic potential for the narrow soliton. 

The freezing of the Bloch oscillations can be easily explained in the adiabatic approximation, where it corresponds 
to the center of the soliton variable being trapped in a local minimum of the effective potential. See Fig. 2| 

Therefore, in the case f > 0, the freezing emerges for solitons centered near lattice sites (and with d^/dt ~ 0); 
while in the case i/ < 0, it emerges when the soliton center is near the middle of an intersite region. In both cases 
the phenomena is favored for narrow solitons (greater w), strong I'-terms (greater v) and weaker external fields 
(smaller |x|). See Fig. [3 The parameter region predicted by the adiabatic approximation [Eqs. 0-®] is in good 
agreement with the results of the numerical evolution of the full equation of motion [Eq. (Q]. However, the adiabatic 
approximation is not able to describe the unfreezing at early times, that can only be observed when the parameters 
are very close to the limit of the freezing region [Fig.|B]. The reason is that this unfreezing is related to non-adiabatic 
effects (change of the soliton shape and radiation) . 



V. CONCLUSIONS 



These results show that inhomogeneity in the lattice is not required to freeze the Bloch oscillations, i.e., to trap the 
soliton in a region with a size of the order of the intersite distance in the lattice. We show that the discrete translational 
invariant integrability-breaking term that makes the transition between AL and DNLS models can freeze the nonlinear 
Bloch oscillations, i.e., gives rise to an intrinsic localization. This phenomenon can be understood as the trapping of 
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FIG. 5: Bloch oscillations arc frozen in the regions above the curves for values of the external field smaller than those indicated 
in the curves (taking for the other parameters 6o = and = 0). Full lines: computation with the complete equations of 
motion. Dashed lines: computation with the adiabatic approximation. 




FIG. 6: Dynamical defrosting due to nonadiabatic effects (radiation and change of soliton shape) of a soliton with initial 
conditions close to the frozen-unfrozen limit in parameter space, wo = 1.2633, = and = in a Schrodinger chain with 
an external field x = 0.3 for v = 0.025. 

the soliton center position in a minimum of the effective potential that we obtain in the collective variable or adiabatic 
approximation. The integrability-breaking terms considered, the ^--term, produces an effective attractive (for u > 0) 
or repulsive (for u < 0) interaction with the lattice sites, that gives rise to local minima in the effective potential where 
the soliton can be trapped. This i/-term has also the effect of changing the main frequency of the Bloch oscillation 
and introducing faster new frequencies in the dynamics, that can also be understood as a consequence of the change 
in the effective potential. 

These new phenomena show a richer dynamics of the DNLS equation, increasing the interest and potential applica- 
tions of this model, and they also suggest the possibility of other new effects in the presence of time varying external 
forces. 
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